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Abstract. The history and recent progresses in the study of bulk viscosity in nuclear 
l/^ I and quark matter are reviewed. The constraints from baryon number conservation and 

electric neutrality in quark matter on particle densities and fluid velocity divergences 
are discussed. 



OO . 

The time scale for damping of the vibration and of the gravitational radiation driven 

^ . instability in neutron stars is crucial to the stellar stability, which is controlled by the 

. shear and bulk viscosities. The bulk viscosity originated from the re-establishment of 

^ chemical equilibrium is important in some circumstances. In this short note, we will 

'p_|. give an overview on the history and recent progresses in the study of bulk viscosity in 

Q . nuclear and quark matter. 

• 1. Nuclear matter 

^> : 

' Sawyer jl] pointed out that at temperatures higher than 10^ K the bulk viscosity for 
^ [ neutron matter is much larger than the shear viscosity. This is determined by opposite 
temperature behavior of the bulk and shear viscosities that the bulk viscosity increases 
while the shear decreases with growing temperature. Since then, the bulk viscosities in 
nuclear matter have been widely investigated. 

npe{fi) matter. Low density neutron star matter is composed of neutrons n and a 
small admixture of protons p, electrons e, and possibly muons /x. The bulk viscosity is 
mainly determined by the modified Urea processes, 

N + n ^ N + P + 1 + Vi, N + p + l^N + n + ui, (1) 

where N denotes the spectator nucleon for energy and momentum conservation, / the 
electron or muon, and i>i the associated neutrino. In late 1960s, Finzi and Wolf analyzed 
the damping of neutron star pulsations via modified Urea processes in npe matter [2], 
followed by Sawyer [1] and Haensel et al. ^ who investigated the bulk viscosity for npe 
and npe/i matter respectively. The relaxation time for the modified Urea processes is of 
order r cx T~^. 
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At densities of a few times normal nuclear matter density, the direct Urea processes, 

n-^p + l + Vi, p + l^n + ui, (2) 

may also be allowed provided the proton fraction exceeds the Urea limit of about 
1/9 [4J. The bulk viscosity was calculated for the direct Urea processes by Haensel 
and Schaeffer [5] in npe matter and by Haensel et al. [6] in npcfi matter. A weaker 
temperature dependence for the relaxation time, r oc T~^, is found due to a smaller 
number of particles involved. Consequently the bulk viscosity for the direct Urea is 
about 4-6 orders of magnitude larger than that from the modified Urea at typical neutron 
star temperatures, T ~ 10^-10^° K. The large difference occurs at low temperatures. 

There may be superfluids in neutron star matter due to attractive part of baryon- 
baryon potential, for reviews, see e.g., Refs. [ZIIH]- The bulk viscosity with superfluidity 
for the direct and modified Urea processes in npen matter was investigated by Haensel 
et al. [31 E], which is substantially suppressed compared to normal state matter. 

Hyperon matter. Hyperons may appear in the neutron star core [9l [IDl [HI fT2] . 
With increasing densities, S~ and A hyperons appear, followed by S~ and S"*". Most 
authors only considered S~ and A hyperons which have the lowest threshold densities. 
Weak nonleptonic hyperon processes such as 

n + n^p + S~, n + p ^ p + A, n + n n + A (3) 

dominate the bulk viscosity, while their direct Urea processes |T3] have negligible 
contributions although they are comparable to nucleons' Urea. The relaxation time 
from processes ([3]) is of order T~^. 

Langer and Cameron introduced the subject of bulk viscosity in hyperonic matter 
by estimating the damping of neutron star vibrations [1A\. Jones carried out the 
semi-quantitative calculation of bulk viscosity in hyperonic matter [T6] . where the 
weak nonleptonic process nn pS^ through VT-exchange was calculated. Recently, 
Jones [17] studied the reaction nn <-> nA, the dominant channel for A production in 
hypernuclei experiments in general, which cannot be mediated by single W boson. 
The rate of the reaction nn nA from the data is found to be several orders of 
magnitude larger than nn ^ pH^ via the VT-exchange. All of these works are the 
order-of-magnitude estimates. 

The various weak nonleptonic hyperon processes have been recalculated by several 
authors with a modern equation of state. Haensel et al. [H] studied the nn <-> pS~ 
process within the nonrelativistic limit, they found the bulk viscosity to be several orders 
of magnitude larger than that of the direct and modified Urea processes. Lindblom and 
Owen [19] computed the contribution of the np ^ pA in addition to the nn ^ pS~ 
process. The superfluidity case had also been studied in both works. 

Most of the weak hadronic processes involved in the nn ^ nA process also 
contribute to the nn one, so the rates of these two processes should be of 

the same order, opposite to Jones' estimate. Dalen and Dieperink [20] studied the bulk 
viscosity from all three processes in ([3]) using one-pion-exchange (OPE), which can well 
describe the rates of these processes in hypernuclei. Their results showed that the bulk 
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viscosity in the OPE picture is about 1-2 orders of magnitude smaller than that with 
ly-exchange. 

Recent hypernuclei data showed the potentials for A and H are attractive but that 
for E is repulsive in normal nuclear matter [211 [22] . Chatterjee and Bandyopadhyay [23] 
took this fact into account and found the disappearance of S in nuclear matter. They 
calculated the bulk viscosity from up pA. Their results showed that the bulk viscosity 
incorporating the hyperon potentials implied by the data is larger than that without 
them. 

2. Quark matter 

When the baryon density is above 5-10 times the normal nuclear matter density, the 
deconfinement phase transition will take place marked by the formation of quark matter. 
Initially there are only u and d quarks, as the density grows the s quarks will appear. 
Quark matter with u, d, and s quarks can be self-bound and be called strange quark 
matter [Ml 125] . The star made of strange quark matter is called quark star or strange 
star [SHE?]. 

Nonleptonic processes. The importance of dissipation due to the nonleptonic 
reaction. 



was first observed by Wang and Lu [28] . They showed that stellar pulsations would be 

strongly damped in quark matter. Saywer [29] then formulated the damping in terms of 
the bulk viscosity based the hnear expansion of the reaction rate of (jlj) in = /i^ — fid, 
where Hi with i = s,d are the quark chemical potentials. In the temperature range 
characteristic of young neutron stars, the bulk viscosity arising from the reaction (jl]) is 
orders of magnitude larger than that for normal nuclear matter. However the above 
linear assumption is not proper at low temperatures (T <^ 6fi), where the rate is 
proportional to 6fi^. Madsen [30| showed that this nonlinearity effect leads to much 
larger the bulk viscosity than previously assumed. Note that strong interactions can 
also influence the rate of the reaction (jl]) and the bulk viscosity of strange quark matter 
significantly [31]. At low temperatures, the viscosity is strongly suppressed, while at 
high temperatures it is slightly enhanced. The above calculations are based on MIT bag 
model where all quark masses are taken to be constants. Alternatively quark massess 
can be assumed to depend on baryon density [32| . By employing this density dependent 
quark model, the bulk viscosity at low temperatures and high relative perturbations 
increases 2-3 orders of magnitude, while at low perturbations the enhancement is 1- 
2 orders of magnitude compared to the results obtained in other approaches [33j . 
Sufficiently dense and cold quark matter is expected to be a color superconductor [31] . 
A recent interest in the study of the bulk viscosity in color superconducting phases is 
growing, for example, the bulk viscosity in the 2-fiavors color superconducting phase 
(2SC) from the reaction (jlj) has been computed by Afford and Schmitt [35]. The bulk 
viscosity due to kaons in color-fiavor-locked (CFL) phase has also been calculated [36] . 




s + u u + d. 



(4) 
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Leptonic processes. Due to phase space restrictions [2ZIEH] the rates of the leptonic 
processes, 

u + e q + Ue, q ^ u + e + 17^, {q = d, s) (5) 

are much smaller than nonleptonic ones at low temperatures. These Urea processes are 
the most efficient way to cool the neutron stars. As shown by Anand et al. [33j, when 
the temperature increases the contribution of leptonic processes to the bulk viscosity 
will exceed that from nonleptonic ones. So more careful studies are needed to explore 
the dependence of bulk viscosity on the broader ranges of temperature. The energy 
emissivity in Urea processes in various color superconducting phases have been studied 
by several groups [39], HOI Hi], |42]. The bulk viscosity from the Urea processes for d 
quarks in a spin-one color superconductor has been calculated by Sa'd, Shovkovy and 
Rischke [13]. 

3. Baryon number conservation and enforced charge neutrality for bulk 
viscosity in quark matter 

A more realistic case for the bulk viscosity in quark matter is to take both of nonleptonic 
and leptonic precesses into account. There must be electrons to compensate the net 
positive charge in a three- flavor quark system because the mass of s quarks is much larger 
than those of light quarks. The bulk viscosity should be calculated for charge neutral 
quark matter. In nonleptonic and leptonic reactions, the electron number and flavors 
are not conserved, while the baryon number and electric charge are conserved. From the 
continuity equations, baryon number conservation and charge neutrality, we can derive 
the following constraints for particle number densities and velocity divergences [44j . 

V ■ = ^^i"^ ■ ' ^eV • Ve = ^ QjUiV ■ Vj , (6) 

i=u,d,s i=u,d,s 

where nj, Qj and are number densities, electric charges and fluid velocities for particle 
j = B,e,u,d, s respectively. The indices B and e denote baryons and electrons. Here 
the first equation is from baryon number conservation and the second one from charge 
neutrality. All previous treatments in the literature used the conditions riBdXi/dt = Ji 
to determine the bulk viscosity, where Xj = ni/ns are fractions of baryon number for 
particle i = u,d,s,e and Jj their sources. Using dns/dt = —ub'V ■ v^, we can get 
UBdXi/dt = dui/dt — Xidus/dt = dui/dt + UiW ■ = Jj. Comparing it with the 
continuity equation for particle i, dui/dt + rijV • Vj = Jj, we obtain V • = V • Vj, 
which obey Eq. ([6]) obviously. Hence the usage of riBdXi/dt = Ji in literature implies 
the unique value of velocity divergences for all particle species. However, considering the 
special role of strange quarks due to their large mass, we propose a new possibility that 
the velocity divergence for s quarks is different from those of light quarks and electrons, 
which corresponds to a new oscillation pattern for the bulk viscosity. As an extreme case, 
we assume V-v^ = and V-v„ = V-v^. From Eq. ([6]), we derive V-Vg = ^^f^r+n^^'^-B 
and V -Vud = ^"^ V-v^. We can use 6nB, Sn^, and as independent variables for the 
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density oscillation. Then 6nu,d and 6P can be expressed in terms of these independent 
variables. In order to close the system of equations for the bulk viscosity, we need two 
additional inputs, e.g. dSrie/dt + ngV • Vg and ddus/dt (note that V • = 0), given by 
the reaction rates which can also be expressed in terms of independent variables. The 
former can be written as dSue/dt + ngV • Ve = ddrie/dt — (2n„ — nd)/{nu + nd)d6nB/dt. 
Then the bulk viscosity is determined from the definition 5P = — CV • v^. The results 
are shown in Fig. 3 of Ref . |13] . The values of the bulk viscosity in this solution is 1-2 



orders of magnitude larger than the conventional ones. 

Note that the assumption V • = reflects the consideration that the s quarks 
are much heavier and may respond to the density oscillation more reluctantly than light 
particles. If the strange quark mass is small, one can investigate many other solutions 
which are close to the conventional solution V- = V- Vj, for example, one can assume 
velocity divergences of particles deviate from that of baryons by a small amount. 
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